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ABSTRACT 


With the aid of the Born approximation, the time 
autocorrelation and power spectral density are calculated for 
the received acoustic signal scattered from velocity fluctuations 
in a turbulent aircraft trailing vortex. The turbulence is 
required to be globally stationary, but only locally homogeneous. 
The treatment includes the effects of spectral broadening due 
to convection of the scattering eddies by a spatially varying 
mean flow and by macroeddies. The 3 dB bandwidth of the 
received signal is related to the scattering angle and the 
core Mach number of the vortex. A primary feature of the 
analysis is that it provides a method for inferring the radial 
intensity distribution of turbulence in a vortex. The analysis 
technique is also applicable to scattering from other turbulent 
flows where significant variations of turbulence level occur 
over distances on the order of the macroeddy size. 
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I. INTRODUCTION 


The detection of aircraft trailing vortices by acoustic scattering 
has in recent years been developed to the status of a practical 
experimental technique. Of notable significance in this area is the 
work reported by Balser, Nagy, and Proudian.^ The pioneering efforts 
of this group have also demonstrated the potential usefulness of acoustic 
scattering as a tool for studying the detailed atvuatiaoe of turbulent 
vortices. The study of vortex structure by acoustic methods is the theme 
to which this paper is addressed. 

The problem to be considered is that of three-dimensional scattering 

from a finite section of a single turbulent vortex. The treatment 

utilizes two theories developed by R. A, Silverman: one his theory of 

local homogeneity, and the other his model for calculating spectral 

broadening due to .macroeddies. The macroeddy model has been applied to 

scattering from globally homogeneous turbulence (as distinguished from 

4 

loaaZty homogeneous turbulence) by Ford and Meecham. While the 
restriction to global homogeneity should be acceptable for describing 
scattering from large volumes of atmospheric turbulence, the resulting 
theories are probably not as useful for predicting signal spectra from 
aircraft vortices and other turbulent flows where the intensity of 
turbulence varies rapidly over distances of a few macroeddies. 

II. MEAN VELOCITY DISTRIBUTION 

It will be evident later on that the correlation time of the received 
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scattered signal is exceedingly small in comparison to observed time 
scales of decay of actual aircraft vortices, which are typically a full 
minute or more. Let us therefore imagine the vortex velocity 
being time-averaged over an interval which is very long in comparison to 
the signal correlation time, but very short in comparison to the vortex 
decay time. Denoting this averaging operation by angular brackets, we 
wish to model the result so obtained by the stun 

( 1 ) 

The spatially varying component represents the mean swirl velocity, 
and Will be discussed presently. The constant vector will represent 
the combined effects of wind and drift in moving the entire vortex as a 
unit. Using (1), we then express the overall velocity field in the form 

where is the fluctuation of ^(^,t) about its mean value. 

The swirl velocity V(r) in (1) will be assumed to be a truncated two- 
dimensional flow having radial symmetry about the z axis in. Fig. 1. It 
will be taken to be entirely circumferential with no radial or axial 
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components. Measurements have showii that the magnitude of increases 
almost linearly from zero at the center to its maximum value, and then 
diminishes more slowly toward zero with radial distance from the vortex 
axis. The maximum velocity attained is commonly called the core velocity, 
which we will denote by V . The radial distance from the axis at which 
this peak velocity occurs, called the core radius, will be denoted by r^. 
The central region of radius r^. is called the core of the vortex. 
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The two-dimensional character of V(r) might be emphasized by writing 
its magnitude as V(a), where a is the perpendicular distance from the z 
axis of Fig. 1. For illustrating the theory, we sLall uje a radial 
distribution of circumferential velocity having a magnitude defined by 

V(a) = f(o) (3) 

where f(0) is a dimensionless distribution given by 

f(a) = Ar^o"V[l - expC- xo® / r^)] (4) 

Definitions (3) and (4) represent a laminar-flow solution obtained 
originally by Lamb,^ and recast here in a fonii similar to that used by 

7 

Georges and others. From the definitions of V and r above, it is seen 

c c 

from (3) that f(cJ) must attain a peak value of tinity, and that this peak 

must occur at a = r . The approximate values of A and x satisfying these 

c 

conditions are A = 1.398 and x 1.256. A plot of f(0) with the argument 

7 

normalized to r is shown in Fig. 1 of Georges' paper, 
c 

It will also be useful to define a vortex core Mach number by 

M = V^/c (5) 

where c is the isentropic sonic velocity in the undisturbed meditun, i.e., 
in the absence of a vortex or turbulence. Even in the largest present- 
day aircraft, the core Mach numbers seldom exceed 0.2. Since the 
macroscale velocity for turbulence in an aircraft vortex should be 
roughly equal to the core velocity, the low Mach numbers involved justify 
a neglect of compressibility effects with respect to all but sound wave 
fluctuations. 
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III. SCATTERING EQUATION 


The scattering configuration is shown in Fig. 1. A uniform plane 
wave is assumed to be incident from below, illuminating a vortex section 
of length L centered about the coordinate origin. The effect of a finite 
antenna beamwidth is taken care of by truncating the vortex. 

Neglecting scattering from temperature fluctuations, assuming 
incompressibility of the flow field, and considering only incident acoustic 
frequencies which are very high in comparison to frequencies associated 
with the turbulence, several authors have correctly derived scattering 
equations from the basic fluid dynamics equations. Of these derivations, 
the om- by Lighthill, which employs the Born (single- scattering) 
app'oicii.iation, is particularly well suited to the present analysis. Using 
Eqs. 6, 7, and 10 from his paper, contracting the tensor indices, and 
introducing a few changes of notation, we have the far- zone expression 

P(^,t) = “ kjx (ip^ c r)"^ (rc)"^ 

X I cos(03 t - k x' - k R*) r*v(r', t - R’/c) d^r* (6) 

In (6), P(r,t) is the deviation from the ambient density due to the 
3 ound wave, replacing the difference (p - p^) used by Lighthill. The 
vectors ^ and replace Lighthill’ s ^ and y, respectively — the unprimed 
vector locating field observation points, and the primed vector locating 
.source points. The vector^ is shown in Fig. 1. The distance R' in (6) 
is the magnitude of the difference vector (r - r’), and the distance r is 
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the magnitude o£ the vector The derivation of (6) incorporated an 
incident cosine plane wave having radian frequency wave number 

- w^/c, and propagating in the positive x direction as indicated by 
the wave vector in Tig. 1. The intensity of this wave, denoted by 
in (6), is given by Eq. 7 of Lighthill's paper. 

The velocity ^ in f6) replaces the velocity fluctuation in 
Lighthill's Eq. 10. Although Lighthill dealt only with turbulent 
fluctuations about the mean flow, it is important to note from his 
derivation that the restriction of v in (6) to representing only 
fluctuations is not necessary for the validity of the expression. 

For the present application, the most appropriate substitution for 
^ in (6) would be the velocity sum in (2) . Retention of the translational 
velocity in (2) would cause a pronounced increase in complication, 
however. Among other things, for sufficiently large values of u , the 
time variation of the scattering angle 0 in Fig. 1 would not be negligible. 
Previous analyses^’® have shown that scattering at different angles arises 
from turbulent eddies of different sizes — which tend to be uncorrelated. 
The signal correlation time from a vortex sweeping by overhead would 
consequently be less than that obtained from a nontranslating vortex; or 
in terms of the power spectral density, a broader spectrum would be 
recorded . 

Even when the velocity u is small enough to neglect the time 
variation of 0 in Fig. 1, the power spectrum would be shifted by an amount 
corresponding to the ordinary Doppler shift from a moving target. This 
latter effect is easily taken into account, however, in an analysis which 
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otherwise neglects Our procedure here will be to substitute (2) into 

(6) with = 0, and then to reinsert at an appropriate point later on. 
Since observed values o£ u in the absence of wind transport. ‘seldom exceed 
1 or 2 ni/sec in magnitude, this procedure should be quite adequate for 
including the vortex drift effect. 

We now must form the lagged mean product ^p(r^,t) p(^, t - T) ^ ^ As 
tlie steps involved in going from (6) to this product axe documented in 
Refs. 4, 8^ and elsewhere, wo simply write down the result and point out 
the minor differences entailed by the present setup. Substituting (2) 
into (6]), neglecting u in accordance with the argument above, and noting 

AA Q 

<uCr.t)) = 0. we construct the product 

pCr, t - T) ) = p^ (2 IT r)"^ c"® cos^S 



where 


K 

= r" 

- r* 

(8) 

rv\ 

yvi 

/wv 

k 

= k 

- k 

i9) 

Aft* 

/viO 


t' 

= t 

- T - R'/c 

(10) 

t" 

= t 

- R"/c 

(11) 


and where R' = Ir - r’ ! and R" = Ir - r’M. The r component V^(r') results 

'/vv /ux ' /lA. I ^ r Va. '' 

from evaluating the dot product in (6) involving the unit vector r/r, and 
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similarly for The scattering vector in (9) is constructed 

from the two wave vectors shown in Fig. 1. These two vectors have equal 
magnitudes, |k^| = |k^| * ® result, has the magnitude 


k = |k| - 2 k^ sin(i6) (12) 

Since the coordinate x in (6) is the x component of the position vector^, 
we note from Fig. 1 that x/r * cos 0. This substitution has also been 
included in (7). 

With the mean vortex flow included, we have obtained two scattering 

integrals instead of one. The first integral in (7) represents scattering 

from the frozen-in refractive index field created by the mean flow, and 

9-11 12 

has been subjected to various two-dimensional and three-dimensional 
treatments. The corresponding power spectral density, defined from the 
Wiener-Kintchine relationship by 


P(r.w) 


(2irPoIi) 


00 

j - -f)) dx 


(13) 


is easily shown to be a line spectrum. 

The sf:Cond integral rises from time-varying turbulent fluctuations; 
its Doppler-broadened spectrum is the subject of the present paper. 


IV. SIGNAL CORRELATION 


We now proceed to evaluate the second integral in (7) . It is 
convenient at this point to regard the velocity component u^(jr^,t) as a 
random process indexed by the continuous parameters (x,y,z,t), and to 
redefine the angular brackets in (7) to mean ensemble averaging. This 
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being done, we assume that the velocity correlation has the following 
form: 

= W(r^) QCJ,4t,r^) C14) 

where C is defined in (8), At is defined from (10) and (11) by 

yvv 

At = t" - t* = T + (R» - R")/c (15) 

and r is the centroidal point 

r. = i(r" + r’) (16) 

It is noted that (14) represents a velocity process which is stationary, 
but not homogeneous (both stationarity and homogeneity, global and local, 
are understood to be wide-sense properties). The dimensionless weighting 
function W(r^) is introduced in preparation for an application of 
Silveirman’s local homogeneity theory; it will renormalize the space-time 
correlation to the local value of mean-square fluctuation velocity. The 
r^ in Q(C»At,r ) arises from our modeling of the frequency spreading due 
to convection of scattering eddies by a spatially varying mean flow and 

by macroeddies. 

8 13 

Lighthill and Meecham and Ford have presented arguments showing 
that the retardation time (R* - R")/c in (15) can be neglected, and we 
accordingly shall replace At in (14) by the signal lag x. Changing 
variables to^and^^, we now denote the second integral in (7) by F(r,T) 
and write 

PCjr.T) = 2irr-‘‘c-'p^I.k;cos^e fw(r,) 

(17) 
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where Re is the real operator, and 

$Ck,T,rJ = (27T)-' |Q(5,T,r 3 e^J^'k (18) 

A^ rM- I A«+ 

In accord with our introduction of local homogeneity assumptions, the 
integral limits in (17) and (18) are extended to infinity in all 
directions.^ 

In order for (14) to be a locally homogeneous covariance as defined 
by Silverman, “ the weighting function W(^^) must be non-negative and 
integrable, and Q(?,T,r^) must be the covariance of a globally homogeneous 
process. Due to the presence of the r^, representing macroeddy and mean- 
flow effects, the function Q does not have the required form, however. 

As noted by Silverman, measurements of Q taken from a system "riding" 
the macroeddies and mean flow would be essentially devoid of large-scale 
influences, and the resulting correlation and spectrum would reflect 
predominantly only the smaller, homogeneous fluctuations. The velocity 
correlation measured from such a system should have the desired local- 
homogeneity form, with the variation of Q with strongly suppressed. 

Although large-scale velocity variations would be essentially 
invisible in such a system, the spatially varying mean-square fluctuation 
level, when normalized to the peak value attained in the vortex, should 
be the same as that measured in the fixed system. We therefore assume 
that W(r^) in (14) is the same for either system, fixed or moving. 

From Ref. 3 we define the velocity of the moving system at each point 
r to be the space-time average of (2) taken throughout a cube centered 
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about r^. The cube must have sides of length much larger than the 
scattering- eddy size (roughly equal to the incident acoustic wavelength), 
but smaller than the macroeddy size. Time averaging is over the 
correlation time of the received scattered signal. We then model the 
average so obtained by the random process 

V (r ) = u + V(r^) + U(rJ (19) 

where u and V(r ) are the same nonrandom velocities as before, and U(r. ) 
is a random process governed by some probability law appropriate to the 
*’smoothed-out" macroeddies. 

3 

Applying Silverman's results to the present situation (see also 

Ref. 4), we have the following relationship between the semitransform 

$(k,T,r. ) in (18), which is measured in the fixed or earth system, and 

its counterpart $ (k,T) in the moving system; 
in 

$(k,T,r^) = {exp[-i k • v„(r J T] V (20) 

where the subscripted angular brackets denote ensemble averaging with 
respect to the probability law of^^(r^) in (19). We next assume that 
the vector components of ^/^^(r^) are jointly Gaussian, with means 
respectively equal to the components of ^^(r^,t)^ in (1), and equal 
standard deviations denoted by U(r^). Performing the ensemble averaging 
in (20) then yields 

exp{-ijcMu^ + V(r^)]T} exp{-4u''(r^^ (21) 

where k is given by (12). We must now detearmine suitable functional 
forms for y(r^) and W(r^). 

>>.+ AV+ 
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While measurements of the relatively young vortices generated in 

14 

wind tunnels have displayed turbulence throughout the entire core, 
observations of more mature vortices in flight tests indicate a laminar 
core.^^ This observed absence of inner-core turbulence is consistent 
with ideas about vortex turbulence advanced by Owen^^ and others. Ov/en 
has postulated the existence of a thin annular region at the outer core 
extremity, i.e., at radius r^, where turbulent energy production occurs. 
The radial distribution of turbulence intensity, and hence our weighting 
S function W(r.), should therefore peak in this region. 

For illustrating the theory, we propose as a qualitative description 
of the radial variation of turbulence level in a mature aircraft vortex, 
the square of f(a) in (4). We then include the effect of finite antenna 
beamwidths in Fig. 1 by truncating the turbulence distribution at 
z = ±L/2, and write in cylindrical coordinates 

W(r^) = g(z3 (22) 

where the gating function g(z) is unity for \z\ < L/2 and zero otherwise. 

Since the macroeddy velocity in a vortex should vary roughly as the 
mean flow, we adopt for the standard deviation >i(r^) the velocity 
magnitude from (3) , 

where g(z) is the same gating function as in (22). 

We now wish to replace $jjj(Jk,T) in (21) by $j^(jk,0) . For a spatially 

4 

uniform macroeddy velocity, an argument given by Ford and Meecham shows 
that for scattering angles not too near 0 = 0 in Fig. 1, the Gaussian 
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factor in (21) diminishes toward zero in x much more rapidly than <I>^(1 ^,t) 
provided that the macroeddy velocity adequately exceeds the velocities of 
eddies whose sizes are on the order of the incident acoustic wavelength 
(i.e., the scattering eddies). In the present situation where the 
macroeddy velocity U(r^.) is allowed to vary in a manner given by (4) and 
(23) , we assume the incident wavelength to be small enough for this 
condition to hold over regions of the vortex which contribute heavily to 
the total scattering — in particular, near the outer-core region. 

Noting from (14) and (18) that 4 is calculated from a correlation 

of r components of velocity, and likewise for in (20), we have from 

17 

Eqs. 3.4.12 and 6.5.3 of Batchelor’s monograph on turbulence, 

*™CJk,0) = k cos^Cie)] / (4iT) (24) 

where a is a dimensionless number equal to about 1.4, e is the dissipation 
rate of turbulent kinetic energy per unit mass, and k is given by (12). 
This expression applies to the inertial subrange of incompressible 
turbulent flows which are also isotropic, and has been cast in terms of 
the scattering angle 0 of Fig. 1. 

For evaluating the integral in (17) , we shall use the cylindrical 
coordinates (a,(^t,z), with measured from the positive x axis of Fig. 1 
in the direction given by the right-hand rule. In this coordinate 
system, we shall denote the component of in the xy plane of Fig. 1 by 
k^, and shall let be the angle between this acj/-plane component and 
the X axis. 

We next regard the mean swirl velocity V(r ) in (19) and (21) as 
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i) 


being in the positive (f direction. With the magnitude of^Qr^) given by 
(3), the definitions of k. and 4 above allow the dot product k»V(r ) 
in (21) to be written 

= - k^V^ f(a) sin«|) - (25) 

We now replace 4 (k,T) in (21) by $ (k,0) in accordance with the 

In /At rn -AAr 

discussion above and substitute (21) into (17). Noting from (24) that 
$^(Jk,0) is a real-valued function, we then use (22), (23), and (25) in 
(17), integrate over z, and write 

F(r,T) = 27rLr‘^c'’®p I. k^ cos^e <b (k,0) 

W •' ^0 10 raW' ' 


exp[-JV^ f^a) k^x^] A(a,x) a da 


(26) 


where r^ is the overall radius of vortex influence, and 

2it 


A(a,x) = J?e{exp[i(w^ - k»u ) x] 


/ 


expEik^Vc 


f(a) X sin(<|> - (j>jj)] 


dtp} 


o 

The integral in (27) has been tabulated (see formula (4.2) of Bowman 
and the result is 


(27) 


A(a,x) = 2TT J^Ek^V^x f(a)] cos E (w^ - x] (28) 

where J is the Bessel function of order zero. 

0 

At this point, it is beneficial to make two changes of variables. 
Since for most functions f(a), the integral in (26) would have to be 
evaluated numerically, it is convenient to normalize a and r^^ to the core 
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radius r^, and to define = a/r^ and second change is 

to replace the time lag x by the dimensionless variable il' defined with the 
aid of (12) by 




Ip = kV^T = 2 k^V^T sin(ie) (29) 

The signal correlatio. s best analyzed in terms of the dimensionless 
ratio F(r,x)/F{r,0) . To construct this ratio, we use the variables 

A\i .Mr 

introduced above, and define 

a 


cm = iN(a^j)] 


■/ 


fV„) exp[-i f'(cr^)i|;^l Jo[i|<k^k 


-1 


“n % 


C30) 


where 

a 

o 

NCV » jf“tv t.„d<T^ (31) 

0 

We then insert (28) into (26), multiply and divide by r^ for normalization 
of cr, and form the ratio F(r,x)/F(^,0) . In terms of definitions (29) 
through (31) , this ratio can be written 

F(r,T)/F(jT,0) s C(iJ>) cost(6o^- Jc‘u^) Tj (32) 

The normalized signal correlation thus separates into the product of an 

envelope function C(ij») and a cosine. 

A plot of C(i/)) for scattering in the xy plane of Fig. 1 (where 

k^ = k) appears in Fig. 2. For this plot, munerical integration in (30) 

was terminated at a = 14, the value of a at which f(cr ) is down to 0.1. 

o n n 

It is noted from (30) that C(ii() is even in \p, and the plot is accordingly 
for positive values of \p only. 


15 



C(\j/} (dimensionless) 


1.0 


0.5 


0 


Figure 2. 



NORMALIZED CORRELATION ENVELOPE FOR SCATTERING IN THE xy PLANE 
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By replacing in (29) by 2 tt/cT, where c is the sonic velocity and 
T the incident wave period, and using the core Mach number defined in (5), 
we can rearrange (29) to give 

t/T = [4ttM sin(i0)] ^ (33) 

From Fig. 2, C(Vi) is down to 0.5 at « 2.3. Vortex core Mach numbers for 
present-day aircraft generally range from about 0.02 to 0.2. For a 
scattering angle 0 of 125 degrees (a local maximum of the scattered 
^ intensity, as we shall see), substitution of these numbers into (33) yields 

a fifty-percent-correlation-level lag of one to ten incident-wave periods, 
depending on vortex size. In order to utilize inertial-subrange eddies 
for scattering, the acoustic frequencies employed would range from about 
1 to 10 kHz. It is thus seen that the signal correlation time is on the 
order of milliseconds at most, which is exceedingly small in comparison to 
time scales of vortex decay for aircraft. 

The sonic velocity c, wave period T, and scattering angle 0 are easily 
determined for most experimental setups i Formula (33) then provides a 
relationship between T, il*, and by way of (5), V^. 

Suppose, for example, the lag at which an experimentally 
determined correlation envelope is down to a given value is measured. A 
plot such as Fig. 2 gives the corresponding value at which the model 
envelope is down to this same value. Substituting these values of and 
into (33), and using (5), then gives the vortex core velocity predicted 
by the model . 

Conversely, given the core velocity of a vortex from which soundings 
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were taken, the value corresponding to a given value could be 
obtained from (33) . It would then be desirable to have several model 
curves such as Fig. 2, obtained, say, for different assumed turbulence 
distributions W^^). A model curve is then sought for which th ‘ value 
C(ifi^) is equal to the experimental envelope value at lag T^, and which 
matches in overall form the experimental envelope. A turbulence 
distribution for the vortex would thereby be inferred. 

These measurements all have analogs in the frequency domain, and we 
proceed now to derive the power spectral density of the received signal. 

V. POWER SPECTRUM 


Recalling that F^>T) is the turbulence contribution to the density 
correlation in (7), we obtain the power spectral density by substituting 
(26) into (13). Before doing so, however, it is convenient to recast 
A(a,x) as given by (27) into a slightly different form. To this end, we 
define 


n = u - k ■ V(r^) - k • u (34) 

o JV\. A\0 ^ ^ 

where k*V(r^) is given by (25). Next, the exponentials in (27) are 
combined, the real part extracted, and, with the aid of (34), the 
resulting integrand is written as cos(J2x). We then use the Euler formula 
to express cos (fix) as a sum of complex exponentials, insert this recast 
version of A(a,x) into (26), then (26) into (13), and write 

^o 2 tt 


P(r,w) 


k^ L (r c)"^ cos^0 


♦ (k,0) 
m'x.’ •' 


I 


f^a) H(<J),o,w) a d()> dcr (35) 


0 o 
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w{ic>re 


HC0,cJ,a)3 


r» 

= i J* exp[-itv2 £*(a) k^T^] dx 


( 36 ) 


In (36), the complex exponentials can be combined, and with the aid of 

to 

formula 861.20 of Dwight, ‘ the integral is 


../j, _ ^ (tt/2) ^ S (o) - £1)^ 

H(<j),a,w) = < exp — r-— r — 

KV^tLo) I [2k^V^fV)J 


exp 


(o) + fl)' 


2 k^ f^(a) 


1 ) 


(37) 


Inserting (37) into (35), it is noted that P(r,w) is even in w. 
Furthermore, since the velocities ^(r^) and appearing in (34) are 
assumed to be low in Mach number, it is noted from (37) that P^,w) in 
(35) consists of two spectral pulses which are centered approximately 
about 0 ) = Since in (37) is likewise a low-Mach-number velocity, 

the two spectral pulses are distinct, having half-power widths which are 
a small fraction of to^. The two spectral pulses therefore make equal and 
essentially nonoverlapping contributions to the total scattered intensity 
as obtained by integrating P(r,t»)) over all values of w. These arguments 
show that it is sufficient to examine only one of the pulses in (37) , and 
we arbitrarily elect the one centered at w “ £1. 

A positive- frequency spectral density could be defined by substituting 
the desired pulse from (37) into (35). Then, by doubling the result, we 
would have a spectral function which would still give the total scattered 
intensity when integrated over u). We accordingly define 

""o 27T 

f 

f(a) exp 


k^ L (27T)'- 


P+(r,w) = 

kV r c 


cos-'e 0 (k, 
in w 


2t 

,0)/ 


(m - £))■ 


L2k^V^ fV) 
c 


. 


o d((i do 


(38) 
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We now introduce three changes of variables in (38) . The first 
change is from (j) to <}»’ = ((j) - ^^). It is noted from (25) and (34) that 
the resulting integrand in (38) would be periodic in <|>', and that the 
limits could therefore be reset to 0 to 2ir as before. 

The second change is from u to r]> where 

n = [w - (w^ - k*u„)]/[2w^M Sin(i0)] (39) 

The variable n is seen to be a normalized deviation from the Doppler- 
shifted incident frequency. It is noted from (5) and (12) that the 
denominator in (39) is simply equal to kV^. 

The third change is to multiply and divide (38) by r^ for 
normalization of the variable a as before. 


We now substitute (25) into (34) , the result into (38) , make the 
variable changes listed above, normalize the result to P (r,0), and 

■t* 


obtain 


a 

0 


J. 

(1 


1 f(a ) Y(n,0 ) a do 

(40) 

P.(T.,0) N,ca„) J 


o 



where 

2ir 

r 



Y(ri.a^) 

= j exp{-4[£(0^)]~^ 

tn - k^k"^ f(o^) sin(}t']^} d4)» 

(41) 


•/ 

0 



and 

( 

'o 

■ 



+ 

/ — \ 

II 

f(0„) y(0,O„) CT^ d0„ 

(42) 


0 


A plot of the normalized spectral density Pjj(n) for the xy plane of 
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Fig. 1 appears in Fig. In this plane, the ratio k^/k in (41) is unity. 
For this plot, the integrations in both (40) and (41) were performed 
numerically, with the limit in (40) set at 14 as done before in (30) . 

It is noted from (40) and (41) that Pj^(ti) is even in n» hence the plot 
for positive values of n only. 

The physical situation is one of small targets — the scattering 
eddies — ■ being blown around by macroeddies and the mean vortex swirl, 
and being translated as a group at velocity^. The various Doppler 
effects due to these motions appear in different parts of (40) and (41) . 

The Doppler shift due to overall translation of the vortex at 
constant velocity is manifested by a shift of the center of the 
spectrum from u»^ to (w^ ” noting that since Pjj(n) 

is even in n, it is centered about n = 0. From (39), then, n “ 0 
corresponds to w = (u> - k»u ). As noted in Section III, the only effect 

* Q At 

of u taken into account in this solution is the spectral shifting. It 

aa-o 

is assvimed that u^ is small enough for its spectral broadening effects to 
be negligible. 

We now wish to obtain a formula for the total 3 dB bandwidth of the 
received signal. Let us denote by n the value of n for which P„(n) = 0.5. 
For this value of n, let the frequency deviation in the numerator of (39) 
be called Ato. Then (39) can be rearranged to give 

= 4Hj.M sin(i0) (43) 

Formula (43) is then an expression for the total 3 dB bandwidth of the 
received signal. 
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(dimensionless) 
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Figui'O 3 





. NORMALIZED POWER SPECTRAL DENSITY FOR SCATTERING IN THE xy PLANE, 
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From Fig. 3, n is about 0.3. For a scattering angle 6 of 125 degrees, 
c 

(43) yields a per unit 3 dB bandwidth of about 1.1 M. In other words, the 
total 3 dB bandwidth of the received scattered signal is expected to be 
roughly equal to the core Mach number of the vortex. 

As a matter of interest, one might inquire as to the relationship 

between Pj,(h) in (40) and C(i/<) in (30) and (32). It can be shown that 
1 . 

[(2'rr)’^ yrid [N^(a^) Pj^(n)l are Fourier transform pairs. 

VI. TOTAL INTENSITY 

Rather than integrating one of the spectral density expressions over 
frequency to obtain the total intensity, it is more convenient to evaluate 
F(r,x) in (26) at x = 0, and to use the formula 

At 

I(r) = c'F(r,0)/p (44) 

Multiplying and dividing (26) by r^ for normalization of o, setting x 
equal to zero, using (12), (24), and (31), and inserting the result into 

(44) , we obtain 

I(r) = I. (rc)‘^ae^'^’7rr“ L N(a^) cos^e cos='(ie) [2 sin(i0)3 

(45) 

This expression has a local maximum at an angle 0 of approximately 125 
degrees, the angle which has been used in numerical examples. 

The angular pattern in (45) and variation as the cube root of 

20 

frequency are familiar from the results of previous investigations. In 
fact, if we were to alter our modeling of the vortex turbulence, and 
regard the illuminated section in Fig. 1 as a "piece" of a globally 
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stationary, globally homogeneous velocity process having a spatially and 
temporally constant mean value, then scattering from the volume of 

4 8 

turbulence so defined would fall within the purview of previous analyses. ’ 

The cylindrical scattering volume, call it is tr^L, or in terms 

of the normalized radial distance a^, is In the globally 

homogeneous case, in (31) is unity, and N(a^) is equal to 

Uien for this special case, the scattering volume can be written as 

V = 2tt N(a ) r^ L (46) 

sc 0 -' c 

Using C46), we can write (45) on a per-unit-volume basis as 

I(r)/V = [I. k* cos^e cotHie) E(k)l / (Sc^r^) (47) 

5 C X O 

where, from Batchelor’s monograph, E(k) = ae ^ k” and k is given by 

O 

(12). Equation (47) is identical to Eq. 25 of Lighthill’s paper, thus 
providing a check on our results. 

VII. CONCLUSIONS 

This paper has presented a model for acoustic scattering from a 
turbulent vortex. The model allows the turbulence intensity, macroeddy 
velocity, and mean flow velocity to vary with spatial location. Although 
the functions used in the present paper to describe these spatial 
variations are in need of experimental confirmation and subsequent 
improvement, they are adequate for demonstrating the method. The technique 
should also be useful for analyzing scattering from other turbulent flows 
where significant variations in turbulence level and mean flow occur over 
distances on the order of the macroeddy size. 
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